Self-diffusion of a sphere in an effective medium of rods 
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Abstract 

Self-diffusion of a sphere in a network of rods is analyzed theoretically. Hydrodynamic interac- 
tions are taken into account according to the model of Dhont et al., under the assumption that 
§ na« 1 and a/L«l, where 1/k is the network hydrodynamic screening length, a = a + D/2, and 

f — 

i ' expressions for the diffusion coefficients are derived and shown to be independent of L. 
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Tracer diffusion of a spherical particle in a suspension of slender particles is a common 
process in various biological systems, such as proteins in suspensions of F-actins, nucleosome 
core particles in DNA dispersions or colloidal particles in polymer systems. 

In Refs. , this process has been used to determine the structure of such net- 

works, in isotropic and nematic phases. A theoretical model of the screened hydrodynamic 
interactions between a tracer sphere and motionless rods has been constructed, by adding 
an effective friction force to the hydrodynamic equations. The model was used to deduce 
from the experiment the screening length of a network of fd viruses as a function of the rod 
concentration. In this communication, we follow the model and derive simple expressions 
for the self-diffusion tensor, helpful in further studies of rod networks. 

According to the fluctuation-dissipation theorem, the self-diffusion tensor T> is related to 
the translational-translational hydrodynamic friction tensor £, 

T> = k B T{Q-\ (1) 

where (...) denotes the average over positions and orientations of the rods. In the above 
equation, the axial symmetry of the system has been taken into account. In this case, the 
translational-rotational hydrodynamic friction tensor vanishes. Therefore the external force 
(F ext ) needed to move the sphere with a velocity v p does not generate a sphere rotation, 
and it is given by the relation, 

(F ext ) = (C) ■ v p . (2) 
In Refs. , the medium velocity v and pressure p are described by the Debye- 



Biiche-Brinkman equation 



A3, 



T](V 2 v - k 2 v) - Vp = 0, V ■ v = 0, (3) 

with the hydrodynamic screening length 1/k, and the isotropic Green tensor in both isotropic 
and nematic phases, 

T s (r) = J-^^I + ^^rr], (4) 
where r = |r|, r = r/r, x = kt and 

hi{x) = -\ + ( 1 + - + \ j exp(-x), (5) 

X 2 \ XX 2 ' 
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h 2 (x) = + l + A )exp(-x). (6) 



0T \ XX 1 



The friction and diffusion tensors for a single sphere of radius a are isotropic, £ = Col, 



Q, 3, 13, 



na << 1, 



T> = P I, with ( = 67rr]a(l + Ka+ (na) /3) |5( and V = k b T /( . In Refs. 
and Co~67r?7a. 

In the network of N rods, the cluster expansion of the friction tensor is performed. After 
the averaging, 

^ = I + ^(M(R,,u,)) + -- (7) 

Co u 

where M(Rj,Uj) is the correction from hydrodynamic interactions between the sphere, 
located at the center of the reference frame, and single rods, labeled j, with their centers 
located at Hj and oriented along iij. 

For the isotropic phase, T> = T>i SO L. For the nematic phase, with the axial direction en, the 
diffusion tensor D = X>||e||e|| +X?_l(I — eyen). With the use of volume fraction <fi = pTcLD 2 /A 
(where p is the number density of rods, D is their thickness and L is their length) and 
Eqs. (CQ), (JTj), the diffusion coefficients V = V iso , T>\\, Vj_ are rewritten as, 

V/V = (l + a^)- 1 , (8) 

with a = ai SO , a\\, a±, respectively, where 

aiso= ^^L J dR ex P[~W R ' u)]TrM(R, u), (9) 

ay is obtained when TrM/3 in Eq. ([9]) is replaced by ey • M ■ en, and the coefficient a± is 
evaluated from the relation 3ai SO = (ay + 2a±). In Eq. (Q, the explicit form of the average 
(...) has been specified, with a direct sphere-rod interaction potential V(R, u), and 

W) = j duf(u)P (u\S), (10) 

where Pq(u|jS') is an axially symmetric probability distribution function, which depends on 
u only through the angle 9, defined by cos 9 = u • ey. Here S denotes the orientational order 
parameter, 

S=p^9- 1 -. (11) 
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To evaluate M in Eq. (J7|), the rods are treated as motionless, and the bead model of each 
rod is assumed, with the odd number of identical beads K — 2n+l — L/D, and their diameter 
D. For a/L« 1 and K» 1, it is useful to represent M as the multiple scattering series, 
with the scattering sequences starting and finishing at the spherical particle; the leading 
contribution to M involves only a single scattering from the rod. This scattering involves 
flows incoming onto each bead —n<i<n, located at rj = R+i/)u, and the flows outgoing 
from any bead —n<j<n, located at Tj = ~R,+jDu, 

n n 

M(R, u)=&^aJ2 • C -(u) • T.fo). (12) 

i=— n j=—n 

In the above equation, the hydrodynamic interactions between the rod beads i and j are 
described by the friction tensors £j-(u), which in general determine the external force Fj on 
a particle j needed to cause the translational velocities Uj of all the particles i, 

F, = ^C ; iu)-f r (13) 

j 

The hydrodynamic interactions between the rod beads are described by the point-particle 
model. Following Ref. jjj, we approximate Eq. f)12p by taking the limit of n — > oo, kD — > 0, 
with L = const, 

-L/2 



(14) 



where R' = R + /u. 

In Refs. [l|, I2I, S] , the hard-core interaction potential is assumed. In this case, 

exp[-/3V(R,u)] = l-x(R,u), (15) 

where x(R> is the characteristic function equal to one inside the cylinder centered at 
R, oriented along u, of length L and radius a = a + D/2, where we used a/L « 1 to 
approximate the shape of the rod close to its ends. 

Eqs. OH]) and (I14p allow to express the hydrodynamic coefficients in terms of four- 
dimensional integrals. In Refs. [l|, Q, , these integrals have been carried out numerically 
for a specific value of T — L/a. In this paper, we will carry out explicitly three of these 
integrations, giving much simpler expressions for the hydrodynamic coefficients. We will 
also derive simple asymptotic expressions in the limit Ka — > 0. 



In the following, we explain our method focusing on the derivation of oti S0 . The other 
hydrodynamic coefficients are obtained by a straightforward generalization. From Eqs. (Q 
and (TT4"j) it follows that we need to calculate 











Tr 


T 

1 s 


(R) ■ (l- 




which is a scalar funct 


ion 


of x = kR' 


and 






A(x, s] 





^(47T?7) 



;A(X, S) 



and 

b (x) = - 4[5hl(x) + Ah x {x)h 2 {x) + 2h 2 (x)% 
b 2 (x)=A[2h 1 (x) + h 2 (x)]h 2 (x). 

For a fixed u, we change the integration variables in Eq. ([9]), R 
Eqs. flUD-CESD, we obtain, 

kq 



ft; 



{KD) 2 \n{nD) Xiso) 



where the hydrodynamic scaling function has the form, 

-3 r L/2 



AL=^r / dl fdR'[l- X (R,n)]A(x,s). 



(16) 

(17) 

(18) 
(19) 

R = R + /u. With 
(20) 

(21) 



The averaging with respect to the rod orientation has been omitted, because in the isotropic 
phase, the integration with respect to R' results in an expression independent of u. 

The hydrodynamic scaling function (12 ip is decomposed into uncorrelated and correlated 
parts, corresponding to the integrands proportional to 1 and (— x), respectively, 



\ h — \ h -l- A 

ll iso — 2 Hso,0 + iv 



h 

iso,c 



(22) 



Calculating A^ o0 , we change R' into spherical integration variables x, s and (f>. We 
obtain, 



dx 



(23) 



The above integral is evaluated analytically by substituting the explicit expressions ( TTBl - 
( [19]) , with hi and h 2 given by Eqs. (JSI)-®, and taking the limit e — > of the exponential 
integrals E n (e), E n {2e). Finally, 



(24) 



The correlated part of the hydro dynamic scaling function follows from Eq. (I2ip . if cylin- 
drical coordinates z', p, <fi are introduced to represent R', with the z' axis along u. In 
addition, we take into account that the integrand does not depend on 0, and therefore, 

.3 rL/2 r l+L/2 



A h 



K' 

2L 



dl 

L/2 Jl-L/2 



I pa 

dz I pdp A{x 1 s). 
Jo 



(25) 



We now change the order of the integration with respect to I and z' in Eq. (125]) . The 
integrand does not depend on I and is symmetric under z'— >— z'. Therefore 

.3 rL 



AL,c = - Y J dz'(L - z') J%dp A(x, s 
and with the change of variables, z' , p —> x, s, 

A h =- 



(26) 



dx Ids + Idx 

JO J KCL J y/l-(Ka/x)'- 



ds 



XS \ 2 , 



1-— )x 2 Ah 
kLJ 



(27) 



The integrand in Eq. ( 1271) is a polynomial in s, therefore the integration with respect to s is 
straightforward and leads to the functions, 



Pn(y) = (n + l) 



ds s r 



i-d-v-) 



2\("+l)/2 



(28) 



Combining Eqs. ff24l) and fl27j) . we express the hydrodynamic scaling function by single 
integrals only, 



20 1 
AL= y + (/i + / 2 )-^(Jo + J 2 ) 



(29) 



where 



In 

J ri 



n + 1 
1 



n + 2 



J dx + J dx p n 
J ' dx + dx p n+1 



(30) 
(31) 



with e = na and n = 0, 2. The above integrals can be evaluated numerically with a high 
precision. 

Assuming a fixed value of kL, we now evaluate asymptotic form of A^ so for e— > 0. Below 
we derive the asymptotic expressions for I ; all the other terms are obtained by the analogical 
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reasoning. To make the presentation more clear, we denote f(x) = bo(x). First, using the 
expansion, 



/(*) = f {0) +f {1) x + o(x), 



(1), 



for x — > 0. 



(32) 



we represent Iq as a sum of two integrals, with fix) replaced by p >+f^x and f(x)—f^ — 
f^'x, respectively. The first integral is carried out explicitly and then the dominant terms 
for small e are kept. Asymptotic behavior of the second integral follows from the expansion 
p n (y) = (n+l)y 2 : /2+o(y 2 ) , valid for y < 1. Combining both integrals, we obtain, 



e-/ (0) + - 
2 J 2 



C U} (kL) 



- In - 

2 J 2 



(33) 



where 



C ot (kL) 



l^-^-^-g + /Oh,«L. 



kL 



(34) 



Asymptotics of the other integrals is also evaluated. I2 has the same form as Iq, but with 
different numerical coefficients of the first and the third terms. J n do not contain a term 
linear in e. Therefore, as expected, in the limit e — > 0, the term (Jo + J2)/{kL) is negligible 
in comparison to (Iq + ^2), which scales as 



A 



h 

iso,c ' 



19tt e 2 
£ — + 2 



32 e 
C {bo+b2} (KL) + — - 161n- 



(35) 



In deriving the above formula, the expansion of the Green functions (jSJ) and ([6]), /ii(x) = 
1/2 — 2a;/3 + o(x) and h2(x) = 1/2 + o(x), has been used to find explicitly the expression 
( |32l) . Then, the function C^^^kL) is approximated by its value at infinity, evaluated 
numerically, 



C {6 ^ 2} (+oo)=-4.944. 



(36) 



For large kL, we analyze the difference, 



C{/}(/tL) - C{/}(+oo) = - / dx 



,r- 



(37) 



With the increasing x, the integrand decays to zero, with terms proportional to 1/x 6 , 
exp(— x)/x 4 , exp(— 2x)/x 2 . For kL > 3, the difference fl37j) is negligible. On the other 
hand, for kL « f, 

For Z^ -* = &q + 62°'' — — 8> an d the above term gives a significantly smaller contribution to 
A^ soc than 197re/4, because e/«L << 1. 

Adding the uncorrected contribution, Af so0 = 20/3, to the correlated one, (1351) with 
(I36p . we obtain a simple approximate expression for the hydrodynamic scaling function in 
the isotropic phase, 

AL = y- — e+[8.41-81ne]e 2 , (39) 

valid for e = Ka < < 1 and d/L << 1 . 

We have applied the same reasoning to evaluate also the hydrodynamic scaling functions 



aiu , introduced in Refs. [l, 2, 3] for the nematic phase. They are expressed in terms of A 
by equations analogical to ( 1201) . with 

A} = A* -25AJ, (40) 
Ai = Af so + ^, (41) 

where 5 is defined in Eq. ( TTTT) and A^ is a function of e = Ka, given by the approximate 
expression, 

Ah U _ 137T + 1 024 2 (42) 

s 15 32 v ; 

The correlated parts of the approximate expressions ( 1391) and ( 1421) (independent of T = 
L/a) are plotted in Fig. [TJ together with the accurate results of numerical integration of 
Eqs. (l29l - (!3~T!) and their analogs for the nematic phase (performed for T = 90.7). For 
net < 0.2, precision of the approximation for A^ is better than 1%, and for A^ better than 
0.3%. 

The specific value of r = 90.7 is cfioseu to match the experiments from Ref. & To 
compare with the results of Ref. [3J], we evaluate the hydrodynamic scaling functions from 
the exact expressions, (1291) and its analog for the nematic phase, and plot the results in Fig. [2, 
together with the numerical results of Ref. [3J|, obtained for S = 0.7 and shown in the left plot 
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FIG. 1: Hydrodynamic scaling functions A^ soc and A^ c versus na. Approximate expressions (|39p 
and (I42p (dashed lines) and results obtained from exact formulas (I29p -(|3ip and their analogs for 
the nematic phase (solid lines). 



of their Fig. 3. The pronounced differences between their and our curves are probably caused 
by inaccuracy of the four-dimensional numerical integration performed in Ref. [3j. In our 
approach, three integrals are carried out analytically, what allows to improve significantly 
precision of the numerical results. In addition, we performed asymptotic expansion of the 
single integral which is left, and proposed simple expressions d3S}-(g2) for A* „ valid in the 
whole range of applicability of the physical model used to describe the system, that is for 
Ka « 1 and a/L << 1. 

Our important conclusion is that in practice, A^ sos depend on Ka only. This result 
indicates that the proper argument of the scaling functions is Ka rather than kL used in 
Ref. Q. 

Possible applications are the following. In Refs. l|, |2|, |3j, diffusion of spheres in isotropic 
and nematic networks of rods has been modelled, based on the isotropic Debye-Bueche- 
Brinkman equation. The model allows to use the measured values of the self-diffusion 
coefficients to determine the hydrodynamic screening length 1/k of the network of rods. 
The simple and accurate expressions for the hydrodynamic coefficients derived in this paper 
can be used to make this procedure efficient. 

In future, predictions of the model from Ref. Q can be verified experimentally by mea- 
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suring self-diffusion of spheres with different radii a in the same network of the rods, and 
checking if k does not depend on a. In the nematic phase, the network of rods should in 
general be described by the non-isotropic Debye-Bueche-Brinkman equation. In the follow- 
ing study, the calculation of the hydrodynamic coefficients described in this paper will be 
modified for such a non-isotropic medium. 
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FIG. 2: Hydrodynamic scaling functions A^, A^ so , A| (top-down). Solid lines: results calculated 
from the exact expressions, ([29]) and its nematic analog; dashed-dot lines: numerical results from 



Ref. 
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